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Abstract 

Recently, Farhi, Goldstone, and Gutmann gave a quantum algorithm for evaluating nand trees that runs 
in time 0(\/N log N) in the Hamiltonian query model. In this note, we point out that their algorithm 
can be converted into an algorithm using 0(N 1 ^ 2+e ) queries in the conventional quantum query model, 
for any fixed e > 0. 

A NAND tree of depth n is a balanced binary tree whose internal vertices represent NAND gates. Placing 
bits x\, . . . ,x%n at the leaves, the root of the NAND tree evaluates to the function f n (%i, ■ ■ ■ ,x 2 «), where 
/„ : {0, l} 2 — > {0, 1} is defined recursively as follows. For n — 0, fo(x) = x, and for n > 0, 

f n (Xl, . . .,X 2 n) = -.(/n-iOl, . . .,X 2 n-l) A /„_i(a; 2 n-i +1 , . . .,x 2 »)). (1) 

The goal of the NAND tree problem is to evaluate f n (%i , • • • , ^2™ ), making as few queries to the bits x\ , . . . , x 2 n 
as possible. The optimal classical randomized algorithm for this problem makes 0(iV - 753 ) queries, where 
N = 2™ [7j [HI [9] . Until now, no better quantum algorithm was known, whereas the best known quantum 
lower bound is only fi(viV) pQ. Here we show that for any fixed e > 0, the quantum query complexity of 
evaluating nand trees is 0(N 1 / 2+t ). 

Very recently, Farhi, Goldstone, and Gutmann proposed a quantum algorithm that evaluates NAND 
trees in time 0(y/N log N), albeit in the unconventional Hamiltonian oracle model [3 [6]. In their version of 
this model, we are given access to a Hamiltonian Ho acting on n + 1 qubits as 

H \b,k) = -x k \^b,k) (2) 

for all b G {0, 1} and k G {0, 1}™, and the goal is to perform the computation using evolution according 
to Ho + Hjj(t) for as short a time as possible, where Hjyit) is an arbitrary driving Hamiltonian (that is 
possibly time-dependent and may act on an extended Hilbert space). 

In the conventional quantum query model, the input is accessible via unitary operations of the form 

U \k,a) = \k,a<Bx k ). (3) 

Two queries of Uo can be used to implement evolution according to Ho for an arbitrary time t, which can 
be seen as follows. The procedure acts on states of the form \b,k,a) (where the last register is an ancilla 
qubit), and consists of the following steps. First, apply Uo to the second and third registers. Then apply a 
controlled- R(t) gate with the first register as the target and the third register as the control, where 

R(t)=( C0S * is[n t). (4) 
V i suit cost J 

Finally, apply Uo to the second and third registers again. With the ancilla qubit initially in the |0) state, the 
net effect of this procedure is the mapping \b, k, 0) i— > cos(xkt)\b, k, 0) +i$m(xf.t)\-ib, k, 0), which corresponds 
to evolution by Ho for time t (that is, the unitary operation e~ lHot ). 
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This simulation of Ho does not imply that any fast algorithm in the Hamiltonian oracle model can be 
turned into an algorithm with small query complexity in the conventional quantum query model. Accurate 
simulation of the evolution according to Ho + Hp (t) apparently requires many interleaved evolutions of Ho 
and Ho{t) each for a small time, yet each of which requires two unitary queries to simulate. Nevertheless, 
it turns out that a Hamiltonian of the kind used in [3] can be simulated in the conventional quantum query 
model with only small overhead. 

In the algorithm of 4 , Hrj(t) is time-independent, so the evolution for time t is given by e - l ( H o+H D )t _ 
Such evolution according to a sum of time-independent Hamiltonians can be simulated using a high-order 
approximation of the exponential of a sum in terms of a product of exponentials of the individual terms. As 
noted in by using a, p th order approximation, the simulation can be performed in 0((ht) 1+1 ^ 2p ) steps, 

where h = max{||iJo||, ||iJD||}, and where the constant depends on p and on the desired simulation error. 
We have \\Ho\\ < 1 and ||-£Td|| < 3, so h < 3. Choosing p to be an arbitrarily large constant, we obtain a 
simulation using 0{t 1+s ) steps, where 5 — l/2p is arbitrarily small. Since the algorithm of [4] applies H for 
time t = 0(\/N\ogN), it follows that the corresponding time evolution can be simulated using 0(N 1 ^ 2+€ ) 
queries to the original oracle, where e > may be arbitrarily small. 

The result can also be deduced by noting that, given query access to the inputs via Uo (Eq. [3]), one can 
easily simulate an oracle for the matrix elements of the underlying Hamiltonian Ho + Hjj used in [1] , and 
then apply results in [2j [3] for simulating sparse Hamiltonians. 
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